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ABSTRACT 

Large galaxies are thought to form hierarchically, from the accretion and disruption of many smaller 
galaxies. Such a scenario should naturally lead to galactic phase-space distributions containing some 
degree of substructure. We examine the errors in mass estimates of galaxies and their dark halos 
made using the projected phase-space distribution of a tracer population (such as a globular cluster 
system or planetary nebulae) due to falsely assuming that the tracers are distributed randomly. The 
level of this uncertainty is assessed by applying a standard mass estimator to samples drawn from 11 
random realizations of galaxy halos containing levels of substructure consistent with current models 
of structure formation. We find that substructure will distort our mass estimates by up to ~20% 
— a negligible error compared to statistical and measurement errors in current derivations of masses 
for our own and other galaxies. However, this represents a fundamental limit to the accuracy of any 
future mass estimates made under the assumption that the tracer population is distributed randomly, 
regardless of the size of the sample or the accuracy of the measurements. 

Subject headings: galaxies: fundamental parameters — galaxies: kinematics and dynamics — galaxies: 
structure 



1. INTRODUCTION 



Measuring the masses of galaxies is fundamental to un- 
derstanding their natu re. B y the mid 1970's, observa- 
in il98(]D clearly indicated that the 



tions (e.g. mubm 
luminous matter in galaxies that we can directly de- 
tect represents just a small fraction of the total mass 
present. This discovery was consistent with theoreti- 
cal work which showed that the stability of the Galac- 
tic disk necessitated t he presence of a previou sly unde- 
tected "massive halo" l|Ostriker fc PeeblesllQTJ ). and be- 
fore long, the now-standard picture of the luminous com- 
ponents of galaxies embedded deep within much more 
extended dark matter distributions was established. 

Today, we have well-developed ideas about the forma- 
tion and evolution of these extended dark matter halos 
— the so-called "hierarchical structure formation" sce- 
nario where small structures in the Universe fo rm first 
and o yer time merge to form l arger structures jPgg^lesI 
ITflfi^l IPress fc SchechteHlTflTJ IBIumenthal et a].in9S4r. 
This picture successfully explains a wide range of ob- 
servations o n larg e scales (e.g. lE iscnstcin et al. 2005; 
iMallcr ct al. -2005'; Tc gmark et"aLl l2004t ISoergel et all 
12003: .Percival et al.. .20021) . but the evolution of lumi- 
nous matter in galaxies within these structures is less 
well understood. For example, in a recent analysis of 10^ 
galaxies taken from the Sloan Digital Sky Survey, the 
mass-to-light ratios of galaxies was f ound to vary sys- 
tema tically with their stellar mass (e.g. lKauffmann et al.l 
1200,'^ . Such trends represent important constraints on 
the general relationship between dark matter halos and 
the luminous galaxies that live within them. 

The trends in dark matter content described above 
were derived from the central properties of a large sam- 
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pie of galaxies. An alternative approach is to use tracer 
populations (such as globular clusters, planetary nebu- 
lae or giant stars) to probe large distances from a cen- 
tral galaxy, allowing much more accurate pictures of the 
mass and extent of the surrounding dark matter halo 
to emerge. This approach derives mass estimates from 
the projected separation and line-of-sight velocities of the 
tracer population. 

A variety of such mass estimators have been developed 
but only a few are commonly applied today. One of the 
earliest is based on the virial theorem and moder n forms 
of it are known as the Virial Mass Estimator l)Zwickvl 
2,933; Limber & Mathews 1960). Its formulation relies 
on independent averages of the kinetic and potential en- 
ergies of the system, resulting in a mass estimate propor- 
tional to (wios^}/(^), where vios refers to the line-of-sight 
velocity of an object, R is its projected distance from 
the center of the system, G is the gravitational constant 
and brackets indicate an average over an entire sample. 
Other estimators, such as the Projected Mass Estimator, 
involve averages over a quanti ty known as the "projected 
mass", defined as (vlM)/G fBahcall fc Tremaineill98ll 
iHeislcr. Trcmaine. & Bahcall 1985). Both these tech- 
niques give estimates for self-gravitating systems and are 
generally applied to galaxy clusters and groups l|Zwickvl 
ll933llT937tlSmthlll9MrLimber fc Mathewslll96fir 

While these techniques have yielded estimates consis- 
tent with the existence of dark matter within galaxies, 
they can only be applied in cases where the number den- 
sity of objects in the system is proportional to the under- 
lying mass density. This assumption is not likely to be 
true for tracer populations around galaxies such as glob- 
ular clusters and planetary nebulae. Evans et al. (200^ 
recently proposed a more general estimator that can be 
applied to these systems, where the densities of the light 
and mass follow different distributions: the Tracer Mass 
Estimator (hereafter TME). 

Irrespective of the application, all of these estimators 
are derived assuming the samples that will be used are 



virialized and isotropic in spatial distribution and ve- 
locity space (or at least described by simple analytic 
functions). At first it may seem that these assumptions 
should be valid for the case of spiral galaxies, since a basic 
model of a spiral galaxy halo is a spherically symmetric 
group of objects on random orbits. Similarly, ellipticals 
may have some degree of flattening but these types of 
anisotropics could conceivably be corrected for when us- 
ing the TME. However, hierarchical structure formation 
suggests that these pictures of galactic structure may not 
be correct in detail. For example, observations of low ve- 
locity dispersions of stars in the around elliptical galaxies 
has recen tly been cited as evidence fo r lack of a dark mat- 
ter halo IjRomanowskv et alJl2003|) . iDekel et. all l|2005D 
show that such low dispersions arise naturally if ellip- 
ticals are formed predominantly from the merger of two 
spiral galaxies (consistent with the hierarchical scenario). 
They found that the orbits of stars at the end of simula- 
tions of the merger of two spiral galaxies tend to be highly 
radial in the outskirts of the merger remnant. This radial 
anisotropy leads to projected velocity dispersions much 
lower than expected at large radii, even in the presence 
of a substantial dark matter halo. 

The hierarchical scenario also naturally leads to struc- 
tures in the Universe forming inside out (i.e. with the in- 
ner parts forming first, see — ^elmi. White, & Springy 
120031 IBuUock fc Johnstoiil 12005(1 . and in sequence in 
scale, smallest to largest. It predicts that structures on 
the scale of galaxy clusters are still forming today, and 
hence that clusters are a poor place to assume virial equi- 
librium, as mass estimators do. In contrast, the main 
epoch of galaxy formation has passed, and the bulk of 
galaxies should be virialized. However, when estimat- 
ing the mass of a galaxy, it is of most interest to find 
a tracer population that extends to large radii to probe 
far out into the dark matter halo. It is precisely these 
regions on galactic scales that are still forming: even 
though the bulk of the galaxy may be virialized, incom- 
plete mixing of debris from recent mergers (i.e. within 
last few Gyr) could lead to significant substructure in 
the phase- space distributions of the tracer popul ation 
' e.g. Johns ton 1998; Hclmi & White 1991 iBullock et all 



120011) . Such substructure would violate the assumptions 
of spatial and velocity isotropy. 

Spatial substructure has been observed in the outskirts 
of our own and other galaxies, lending strong support for 
the hierarchical structure formation model. For example, 
when Ibata et al. (1994, 1995) discovered the Sagittarius 
dwarf galaxy it was clear that it is currently merging with 
the Milky Way. While this work only hinted at "tidal dis- 
tortions" , later work by various groups found associated 
streams of debris ext ending all around the Galaxy (see 
iMaiewski et al.ll2003l for a summary). The tidal debris 
of other satellites of the Galaxy h ave also been inapped 
(e.g. the Carina dwarf spheroidal 'Ma ie wski et ani200ff l 
and similar features hav e been di scovered in the halos of 
other galaxies (e.g. M31 llbataet alJI200H) . 

As noted above, samples with spatial overden- 
sities or correlated motions would break the as- 
sumption of a tracer population randomly dis- 
tributed in phase-space. For example, as globu- 
lar clusters and gi ant stars can be used to dis- 
cover tidal features llLvnden-Bell fcLvnden-Bell 1 119951 
iPalma. Maiewski. fc JohnstonI I2002D . it is clear that 



they may not be distributed randomly, even in samples 
where spatial subst r ucture is not immediately obvious. 
iWilkinson fc Evana l|1999j) estimated quite how badly 
these broken assumptions might affect mass estimates by 
comparing the result of applying a mass estimator to a 
truly random test data set with one in which all the data 
lies along two streams. They concluded that there will 
be a 20-50% error in the mass estimate in this particular 
case, considering it an upper limit since the substructure 
example they chose was so far from random. However, 
this work did not look in general at the expected level of 
substructure around galaxies today. 

The goal of this paper is to assess the effect of sub- 
structure on mass estimates for galaxies formed hier- 
archically in a consistent cosmological context. Rather 
than looking at spe cific, contrived examples of substruc- 
ture (such as Wilki nson fc Evansll 999'l we will draw sam- 
ples from m odel stell ar ha. los built entirely by accretion 
i|Bullock fc Johnstonl l2005|) and use the TME to provide 
mass estimates for these halos. These samples were cho- 
sen to be similar to true observations of the different 
types of tracer populations, allowing us to assess the ac- 
curacy and precision of estimates that have been made 
and others that could be made in the future. 

The TME and model datasets are discussed in greater 
detail in ^ while the samples taken from these datasets 
and the mass estimates provided from them are described 
in ^ Our conclusions are summarized in §4. 

2. METHODS 

2.1. The Tracer Mass Estimator 

To derive their estimator. lEvans et al.l l(2flfl,'Jl assumed 
that the tracer populations and underlying mass distri- 
bution of a system can be approximated by certain stan- 
dard forms. The tracer population itself, gathered be- 
tween some spherical radii Tin and rout , must have a num- 
ber density that is spherically symmetric and described 
by a single power law in this region: 
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where po is the density at r = a or at all radii when 
7 = 0. The underlying potential of the system is also 
characterized by a general form: 
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where 6 is a constant with units of length and vq is the 
constant circular velocity of test particles when a = 0. 
The form of the potential for a = describes galax- 
ies with flat rotation curves, for a — 1 the potential is 
Keplerian, and a — j — 2 describes a self-gravitating, 
mass-foUows-light system. 

Assuming a tracer population, with density distribu- 
tion given by equation (1) and moving with an isotropic 
velocity dist ribution in a poten tial of the form given in 
equation f2 HEvans et alJ 1(20031 ) determined how the av- 
erage value of the projected mass (vf^^R) relates to the 
mass within the outermost radius 



M(ro 



\ ^ / 2 



N 



C 1 



i—1 



los,i 



Rj 



(3) 



where N is the total number of objects in the sample and 
C is defined as 
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(4) 

These last two equations define the TME. Note that 
while 1^ and ^ assumed the existence of the unknown 
constants po, vq, a, and b in order to derive the above 
equations, no knowledge of the true value of these con- 
stants is necessary in order to use this estimator. 

Evans ct al. (2003) also derived the TME for systems 
with anisotropic velocity distributions in which the ra- 
tio of the radial to the tangential velocity dispersion is 
allowed to vary. We consider only the isotropic version 
since we expect the uncertainty due to substructure to be 
similar with either form of the estimator. Substructure 
represents a deviation from a random distribution that 
cannot be described by simply introducing anisotropy to 
the velocity distirbution. 

2.2. Discussion of Parameters 

The primary advantage of the TME is that it can be 
adjusted for different populations, giving it the widest 
possible application. This is apparent in the four differ- 
ent parameters found in equation Q for C: a, 7, rin, 
and Tout- Since C is the proportionality constant for the 
mass estimated from a given set of observational data, 
determining proper values for these parameters is vital 
to retrieving accurate estimates. 

The parameter 7 is the power law fall-off of the den- 
sity of the tracer population, and can be estimated di- 
rectly from the tracer's surface density distribution in 
the region of interest (which should be chosen so that the 
tracer population is well-represented by a single value of 
7). Tests deriving 7 for different random samples drawn 
from the same population found variations in the esti- 
mates of order 12% for a sample size of A^ = 50, 3% for 
N = 500 and <1% for N = 5000. 

The radii rjn and rout represent the inner and outer 
spherical radii of the sample, and hence are unobservable. 
However, since extended samples over a large range of 
radii are of most interest, and the density of tracer pop- 
ulations typically falls rapidly with radius (e.g. 7 ~ 3 
for the Milky Way globular cluster system), the spheri- 
cal radii should be reasonably approximated by the pro- 
jected inner i?in and outer -Rout radii of the sample. 

The most uncertain parameter is a, which represents 
the unknown shape of the underlying potential (in fact, 
if the potential of the system were somehow known, then 
using a mass estimator would be irrelevant). It is true 
that the behavior of the potential only needs to be known 
in the region prob ed by the tracer population, which 
IE vans et al.l (|2003f) always intended to be the outer re- 
gions of galaxies, where values of a « (corresponding 
to a flat rotation curve) should be appropriate. (Indeed, 
l^ilkinson & Evans 1999, found that this approximation 
is valid out to nearly 170 kpc in the Milky Way.) How- 
ever, A^-body simulations of structure formation suggest 
that dark matter halos in general should have potentials 
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Fig. 1. — The values of a corresponding to different values of 
— for NFW halos, where Tb is the scale radius. 
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IjNavarro. Frenk. fc Whit^ 1199^ iBullock et all l200l(l . 
where Mhaio and Tg are the mass and length scale of the 
halo respectively. Figurenshows what values of a are ex- 
pected for various regions of an NFW potential. For very 
large/small radii, the values of a tend toward +1/ — 1. 
In practice, Klyoin ct al. (2002) estimate rg in the range 
15-30 kpc in order to fit observations of the Milky Way 
and Andromeda, so our interest lies in the middle region, 
where r is on the order of 0.1-10 rg (i.e. few to hundreds 
of kpc). No single value of a would seem appropriate 
here and, in fact, in the region around r~rsa = — 0.3 
rather than a = 0. 

In summary, the values of rju, and rout are known with 
some confidence because they can be estimated from the 
tracer sample itself. The uncertainties in 7 and a pro- 
vide the main source of error in C (and consequently the 
mass estimates). Figure El shows C as a function of a 
for three different values of 7. The greatest range in C 
occurs for the highest value of 7 — tracers whose den- 
sity falls much more rapidly than the underlying mass 
distribution — where large adjustments (i.e. large val- 
ues of C) to the raw mass estimates are necessary. For 
an observed sample we might expect 7 ~ 3 (e.g., the 
globular cluster system of the Milky Way) and a to vary 
in the range —0.5-0.0 in the region of interest (as dis- 
cussed above and illustrated in Figure 0. The range in 
a leads to an uncertainty in C (and the mass estimate) 
of order 10-20%. For small sample sizes (A^ < 50), we 
showed that we will be unsure of our estimate of 7 at 
the 12% level, which adds another ~ 10% uncertainty to 
C. For larger samples, 7 will be better known and this 
contribution to the error budget will become negligible. 

2.3. Model data 

2.3.1. Description 

Our mock tracer population samples are drawn 
from the eleven s t ellar halo models simulated by 
IBullock fc JohnstonI l)2005f) (hereafter referred to as Ha- 
los 1-11). These models were created within the context 
of a standard A Cold Dark Matter universe, and hence 
are expected to contain a degree of substructure rep- 
resentative of true galactic halos. IBullock fc JohnstonI 
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Fig. 2. — C as a function of a with Hn/rout = 20 for three 
different values of 7 



(J2005J) began with a Milky Way potential consisting of 
three components: a halo, a disk, and a bulge. Each of 
these components grows over time according to analytic 
formulae. More than a hundred N-body simulations of 
dwarf galaxies disrupting within these potentials are run 
for each halo, with the unbound particles contributing to 
the growth of the halo. The eleven models result from 
eleven different randomly-realized accretion histories. It 
is these histories that dictate the growth of the parent 
galaxy and the accretion time and mass of each satellite. 

Note that the stellar halo is represented in these models 
by assigning a variable mass-to-light ratio to each dark 
matter particle in the simulations. These mass-to- light 
ratios are chosen so that the luminous matter in the in- 
falling dwarfs initially follows a King model embedded 
deep within a dark matter halo in the shape of an NFW 
profile. The samples described in ii3.1l are drawn at ran- 
dom from the full set of particles, weighted inversely by 
their mass-to-light rat i os. 

B ullock fc JohnstonI l)2005f) confirmed the consistency 
of their models with the size and density profile of the 
Milky Way's stellar halo and the number and distribu- 
tion in structural parameters of the Milky Way's satel- 
lite population. None of the resulting halos are perfectly 
smooth; they all contain some amount of substructure, 
which varies due to the different accretion histories. In 
addition, each halo is somewhat flattened in response to 
the disk component of the parent galaxy potential. 

2.3.2. Isolating effects of substructure 

Although we are specifically interested in the effects of 
substructure, our models suffer from additional limita- 
tions that compromise the accuracy of the TME. First, 
the flattening described in H2.3.1I is a large scale devia- 
tion from isotropy that will clearly affect mass estimates: 
if the datasets are viewed "face-on" {i = 0°) with respect 
to this flattening the velocity dispersion (and hence mass 
estimates) will be smaller than if they were viewed "edge- 
on" {i = 90°). Second, because the models are created 
by superposing individual simulations of satellite accre- 
tion rather than running a single, fully self-consistent 
simulation, we cannot expect them to be virialized. 

The incorrect assumptions of spherical symmetry and 
virialization mean that we cannot expect the TME to 
recover the true mass of the underlying parent potential 
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Fig. 3. — The observable properties and projected positions 
of a 5000 particle sample taken from Halo 5 with an inclination of 
90°. The top row shows an azimuthal angle of 0° while the bottom 
shows 90° . Differences in the location of substructure in the Vi^^ 
vs R panels suggest differences in mass estimates made using this 
data set. 



in our simulated data sets, even if substructure were not 
present. Hence, rather than comparing to the true mass, 
we isolate the effects of substructure on our estimates 
for a given halo model by looking at how they vary along 
different lines of sight (i.e. azimuthal angle) all chosen at 
a fixed inclination to the plane of the parent galaxy disk. 
In a system with overdensities and streams, the corre- 
lated features would appear to move around the system 
(in R and wios) when different angles are probed, as il- 
lustrated in Figure |21 which will show up directly as a 
varying mass estimate. 

3. RESULTS 

3.1. Simulated samples 

Samples are drawn from the test data set in two dif- 
ferent ways to simulate: (i) an extra-galactic globular 
cluster survey (i. e. that might be conducted today; 
IE vans et al.1 1200!^ : and (u) a large data set of giant 
stars (that might be collected in the future around Local 
Group galaxies). 

3.1.1. Simulating globular cluster surveys 

For this example, we assume that globular cluster 
systems are formed primarily through the accretion of 
satellite galaxies and their associated globulars. This 
type of chaotic forma tion scenario has been proposed by 
ISearle fc ZinnI ljl978|) . who suggested that globular clus- 
ters in the Galaxy's outer halo originated in infalling pro- 
togalactic fragments. Indeed, observations indicate that 
at least some of the Milky Way globulars were likely to 
have been accreted from smaller galaxies, since the Sagit- 
tarius dwarf is apparently con tributing four glo bulars 
to the overall Galactic system l|Ibata et aLlllQQ^^ . Fur- 
thermore, the properties of the globular clusters in the 
Galaxy's outer halo — such as their horizontal branch 
morphologies, kinematics, and overall spatial distribu- 
tion — seem to suggest that the y were formed in satel- 
lite systems that were accreted IjZinn 1 11993') . On the 
other hand, a number of other scenarios have been pro- 
posed for how globular cluster systems form (e.g., Forbes, 
Brodie, fc G rillmair 1997, Co te, Marzke, & West 1998). 
For example. I Ashman fc Zepa (1992) theorized that a sig- 
nificant fraction of a galaxy's globular clusters could be 



formed from dissipating gas during galaxy-galaxy merg- 
ers, and observations of nearby merger events indicate 
that massive star clusters do form in this way (Whitmore 
2000). Therefore the results in this section are likely to 
represent an upper limit on the effect of substructure on 
mass estimates using globular clusters as the tracer pop- 
ulation. 

Selecting particles to represent globulars 

In order to select a realistic sample of N globular clus- 
ters from our models an initial cutoff is made such that all 
satellite galaxies with total luminosities less than IO^Lq 
are removed from consideration, as it would be very 
unlikely that they would contain any globular clusters, 
which themselves may have luminosities from W^Lq to 
~ 10^ Lq. For the remaining satellites we have made 
the simplifying assumption that the number of globu- 
lar clusters is proportional to the satellite's luminosity 
— i.e., the specific frequency of globulars is constant. 
In fact, the relationship between the properties of dwarf 
galaxies and the number of globular clusters they host is 
not well-determined, and (for example) the Local Group 
dwarf galaxies exhibit a wide range of properties in terms 
of the relative number and metallicity range of their glob- 
ulars llGrebe]l B)021. 

Next, given a line-of-sight to "view" the model and the 
range of projected separations of the globulars from the 
parent (see below) all particles from satellites above our 
threshold luminosity value within the specified range of 
separations are selected and saved. Each satellite is then 
assigned a sample size, which is a fraction of the desired 
size N, in accordance with our specific frequency assump- 
tion. These samples are gathered from the individual 
satellites by designating particles at random (weighted 
by the luminosity associated with each particle — see 
§2.3.1) to be globular clusters. Combining these gives 
the final sample. 

To illustrate the results of the selection process. Fig- 
ure^ shows how the luminosity, and thus the globulars, 
will be distributed among all eligible satellites in the 
5 < R < 100 kpc regions of two different halos chosen to 
be at the extremes of the distribution of properties for 
all eleven halos. In these examples. Halo 1 has 16 satel- 
lites accreted ^ 8-10 Gyr ago contributing to this region 
with the most luminous satellite containing about 16% 
of the globulars. In contrast. Halo 9 would have only 3 
contributing satellites, 2 of which were accreted 4-5 Gyr 
ago, with the most luminous containing nearly 75% of 
the globulars. 

Given the extremes apparent in FigureQl Figure[Slsum- 
marizes the general properties of infalling satellites for 
all halos by showing the distributions for the number of 
eligible satellites and largest fractional luminosity con- 
tributed to a halo by any single satellite. The variety of 
properties for the infalling satellite population evident in 
Figures 01 and is due to the different accretion histories 
of each halo model. Overall, we expect the level of sub- 
structure in a halo (and hence a globular cluster system) 
to depend on how many satellites contribute significantly 
to the halo's mass, how massive the largest contributer 
is, and how far ba ck the most recent merger took place 
(see [Bullock & ,Toh nstonll2005l for a discussion). These 
figures suggest that the TME may perform differently for 
the different models. 
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Fig. 4. — The fraction of globulars assigned to the eligible 
satellites (i.e. satellites with luminosities greater than IO^Lq) for 
two different halos, plotted as a function of the satellite's accretion 
time. 
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Fig. 5. — The distribution of the number of satellites used (upper 
panel) , and the largest fraction of globulars contributed by any one 
satellite (lower panel) for all eleven halo models. 



Number and extent of globular cluster samples in spiral 
and elliptical galaxies 

The range of radii around the parent galaxy over which 
the sample will be taken and the number of globulars that 
you would expect to be found in this range depends on 
the morphological type of the parent. 

Nearly half of the Milky Way's ~ 180 globulars lie 
within 5 kpc, while the mo st distant is located at around 
100 kpc IjAshman fc ZeDil[l998j) . In other spiral galax- 
ies, most are located within around 20 kpc from the cen- 
ters of the host galaxies (Rhode, private communication). 
Given these ranges, the systems in spirals are represented 
in our datasets by samples of 50 globulars taken from 5 
kpc to 50 kpc. 

Elliptical galaxies typically have more populous glob- 
ular cluster systems than spiral galaxies. Some giant 
ellipticals have globular cluster systems with thousands 
of clusters, extending well beyond the light of the galaxy 
(e.g., Ashman & Zepf 1998, Rhode & Zepf 2001, 2004). 
For our simulations, we chose 500 particles over the ra- 
dial range of 5 to 100 kpc. It should be noted here that 
the halo models are specifically for spiral galaxies, not 
ellipticals; ellipticals are likely to have undergone much 
more violent, recent mergers which could effectively mix 
any phase-space substructure resulting from minor merg- 
ers. Thus, applying the results of this study to ellipticals 
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will give us an upper limit to the size of the error in the 
mass estimates as it will overestimate the level of sub- 
structure existing in these galaxies. This example also 
serves to give us an indication of how the uncertainty in 
the mass estimates due to substructure is related to the 
sample size of a tracer population. 

3.1.2. Simulating giant star surveys 

Using systems of giant stars in extragalactic halos of- 
fers an alternative to using globular clusters. Such stars 
can be selected photom etricall y with some confi dence 
l|Maiewski. Ostheimer~K unkcl. fc PattersonI l200f]D . Al- 
though the number of spectra are currently limited to a 
few hundred stars in fields probing a small fraction of the 
halos of Local Group galaxies, it is feasible for such sam- 
ples to be increased dramatically in both size and area 
coverage in the near future. 

In this case, there is no luminosity cutoff that needs to 
be made and no assumptions concerning the frequency 
with which these objects would be found in the satellite 
galaxies. The particles are chosen from a range of 5 kpc 
to 100 kpc, which allows the resulting estimates to be 
compared to those from the samples described above. A 
number of 5000 is chosen as the sample size. While this 
is significantly larger than sample sizes of surveys con- 
ducted today, it is also much less than the total number 
of giant stars that exist in the halos of spiral galaxies. 

3.2. Mass estimates 
3.2.1. Uncertainties due to substructure 

Each panel of FigureElshows the results of applying the 
TME to samples drawn from Halo 1 at fixed inclination 
(indicated in the panel) and as a function of azimuthal 
angle. (In this, and all subsequent figures, Mttuc is the 
known total mass in the simulation — galaxy plus dark 
matter halo — enclosed within the outermost projected 
radius of the data sample.) A value of a = was as- 
sumed for every sample, while 7 was calculated from the 
surface density of the data. Each point represents the 
mean value of 100 different samples at the given inclina- 
tion and azimuthal angle, while the error bars give one 
standard deviation. As discussed in ^2.'S.2\ we do not ex- 
pect the average estimate to give a correct value even in 
the absence of substructure: the datasets are not virial- 
ized or spherically symmetric and a is unknown. Rather, 
we interpret the scatter in the mean values in each panel 
as a measure of the effect of substructure. 

The top row of panels in Figure shows the mass esti- 
mates for our simulated globular cluster survey of a spiral 
(the most relevant sample to current observations). The 
scatter in the mean lies well within the error bars for all 
points in this row (i.e. across both azimuthal and incli- 
nation angle). Hence the uncertainty in the estimate is 
clearly dominated by small number statistics. Moreover, 
our analysis excludes individual uncertainties in the ob- 
served positions and velocities, which have been shown 
to have a sign ificant impact on th e final uncertainty of 
the estimates (| Wilkinson fc Evaiiall9 99'l. These results 
generally hold for our full sample of 11 halos, although 
some halos have differences between values that are just 
beginning to reach a notable level of significance when 
comparing between inclinations of 0° and 90°. Overall, 
it is safe to say that for sample sizes associated with 



current observations, substructure does not contribute 
significantly to the error budget of mass estimates made 
with the TME. 

The lower two panels of Figure|Sldemonstrate that vari- 
ance in mass estimates due to substructure effects should 
become a dominant contributor to the error budget once 
sample sizes of order 500 or greater are obtainable. It 
is interesting to note that the points (i.e. mean values 
of the mass estimates) in the middle and bottom panels 
fall at similar locations, which suggests similar results 
from our elliptical globular cluster and giant star samples 
even though the former is restricted to being drawn only 
from debris from the more luminous satellites. In fact, 
^80% of the mass in these stellar halos comes from th e 15 
most luminous satellites (jBuUock fc Johnstonll2005F l. so 
the different sampling methods actually produce similar 
samples. 

Figure[3summarizes the results for all halos — the his- 
togram shows the distribution of the standard deviation 
of mass estimates made from samples of size N = 5000 
taken at 51 equally spaced azimuthal angles at an in- 
clination of 90° for each halo. Based on this analysis, 
the effect of substructure can be quantified as introduc- 
ing an uncertainty in mass estimates in the range 0-20%, 
depending on the accretion history of the galaxy in ques- 
tion. This represents a fundamental limit to the accuracy 
and precision of measuring the mass of a galaxy using the 
TME which can not be overcome by increasing the sam- 
ple size. 

3.2.2. Uncertainties due to oblateness or triaxiality 

Using the larger samples, it also becomes apparent that 
the intrinsic oblateness of the stellar halos causes the 
mass estimates to vary along different lines-of-sight. This 
is demonstrated in Figure |Sl which shows the mass esti- 
mate as a function of inclination of the line of sight to the 
disk for our most flattened stellar halo, Halo 3 (minor- 
to-major axis ratio c/a ~ 0.6). The points in each panel 
vary little in azimuthal angle (i.e. due to substructure), 
but the average mass estimated at an inclination of 90° 
is nearly twice that of the average mass given for an in- 
clination of 0°. 

If we can assume that the stellar halo of a spiral galaxy 
is purely oblate and that its symmetry plane is aligned 
with the disk, the systematic uncertainty in the mass 
estimate due to inclination to the line-of-sight could be 
accounted for by measuring the apparent flattening of 
the tracer population and correcting by the known in- 
clination of the disk. These assumptions arc not too 
unreasonable. Although dark matter halos that form 
in cosmological N-body simulations tend to be triax- 
ial, with minor-to- r najor axis ratios of order cla ~ 0. 6 
(l.ling fc Sutcll200l lB7aiocll2nol lAflgood et al.ll200l . 
simulations with a baryonic component, including the 
effects of gas cooling and star formation lead to only 
mildly triaxial halos ic/a y 0.8 and b/a ~ 0.9, see 
iKazantzidis et a l.''2004l' 2005HBailin et alJ2005f) contain- 
ing disks whose rotation ax es align with the halo's mi- 
nor axis l|Bailin et alJl2005^ . Hence, after correcting by 
the disk inclination, the remaining asymmetry should be 
very mild {b/a ~ 0.9). Our own results (i.e. factor two 
variance in the mass estimate for c/a ^ 0.6) suggest that 
this would lead to an uncertainty in mass estimates of 
order 10%. 



4. CONCLUSIONS 

This paper presented an analysis of errors inherent to 
the Tracer Mass Estimator, the most appropriate estima- 
tor to use when calculating the mass of a single galaxy 
using tracer population kinematics. It was applied to 
samples drawn from 11 halo models which contained lev- 
els of substructure consistent with the current A-CDM 
paradigm of galaxy formation. Each of the following were 
found to contribute to the error budget for the mass esti- 
mates at the 10% level: (i) incorrectly assuming a scale- 
free power law potential throughout the region sampled; 
(ii) uncertainty in the measurement of the surface density 
distribution of the tracer population; (iii) oblateness or 
triaxiality of the tracer population; and (iv) substructure 
in the tracer population. (In our eleven examples, sub- 
structure caused misestimates in masses of up to ~20%.) 

This level of systematic uncertainties is similar to 
that expected for mass estimates of galaxy clusters us- 
ing intergalactic planetary nebulae. In her anal ysis of 
a num erical model of a galaxy cluster, Willma n et al.l 
^^2004^ found that applying the Projected Mass Estima- 
tor (PME) to unbound particles at the end of the simula- 
tion led to an overestimate of the cluster mass by of order 
10%. Her study suggests that this uncertainty was due 
to a combination of anisotropics in the velocity distri- 
bution of the tracer population with known limitations 
of the PME, rather than substructure. However, since 
the study was limited to one galaxy cluster, it could not 
address the effect of substructure more generally. 

Systematic errors of up to '^20% are not the main con- 
tributors to the error budget when using smaller sam- 
ple sizes consistent with current surveys of spiral galaxy 
globular cluster systems {N < 100). In these cases the 



statistical and measurement errors dominate. However, 
once sample sizes of order A^ = 500 are achieved (either 
galaxies with larger globular cluster systems or future 
samples of giant stars in nearby galactic halos), these 
systematic errors represent a fundamental limit of the 
TME. 

Of course, a larger sample size can be utilized to re- 
fine the TME to solve for additional global properties, 
such as the full density profile of the underlying mass 
distribution (rather than assuming a single power law) 
and simple asymmetries in the tracer population distri- 
bution (i.e. to correct systematic errors (i), (ii) and (iii) 
above). Moreover, in the case of the Milky Way, future 
astrometric missions such as NASA's Space Interferom- 
etry Mission and ESA's Global Astrometric Interferom- 
eter for Astrophysics promise additional dimensions of 
phase-space information at very large distance from the 
Galactic center so that phase-space anisotropics could 
be directly assessed. However, substructure in the tracer 
population is not something that can be simply parame- 
terized and hence our study suggests a fundamental limit 
of 20% to the accuracy of any mass estimator that is de- 
veloped under the assumption of randomness, however 
large the sample or accurate and complete the measure- 
ments. 
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Fig. 6. — Mass estimates for Halo 1 as a function of aziniuthal angle and at the specified inclination of the line of sight . (Here, i = 0° 
is face-on and i = 90° is edge on.) The top row simulates a globular cluster survey of a spiral galaxy with A^ = 50 and 5 < -R < 50 kpc. 
The middle row simulates a globular cluster survey of an elliptical galaxy with N = 500 and 5 < R < 100 kpc. The bottom row simulates 
a giant star survey with N = 5000 and 5 < i? < 100 kpc. 
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Fig. 7. — Histogram of the distribution of uncertainties due to substructure for all eleven halo models . The mean of the quantity 
Mestimate/^true was Calculated from 10 A'^ = 5000 samples at each of 51 azimuthal angles equally spaced in the range 0° to 180°. The 
uncertainty was quantified by taking the standard deviation of the 51 means. 
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Fig. 8. — Same as the bottom row of panels in Figure|^but for Halo 3. The eff'ect of inclination is seen in this case as a clearly increasing 
mass estimate. 



